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ABSTRACT 


A "constraint length" for a convolutional code to be used with sequential 
decoding is usually specified by the allocation of a certain number of bits 
to be used in the "tail" of an encoded frame. The constraint length, K, of 
the code has conventionally been chosen to match this tail parameter; this 
report shows that several other options are available. Codes in various 
options which require the same tail allocation are said to have the same 
"effective constraint length", K^, regardless of their actual constraint 
lengths. 

Virtually all previously-suggested rate 1/2 binary convolutional codes 
with = 24 (which is the effective constraint length specified for the 
International Ultraviolet Explorer spacecraft) are compared in this report. 
Their distance properties are given; ar.d their performance, both in com- 
putation and in error probability, with sequential decoding on the deep- 
space channel is determined by simulation. Recommendations are made both 
for the choic of a specific = 24 code as well as for codes to be included 
in future coding standards for the deep-space channel. 

A new result given in this report is a method for determining the sta- 
tistical significance of error probability data when the error probability 
is so small that it is not feasible to perform enough decoding simulations 
t"* obtain more than a very small number of decoding errors. This result 
should be of general usefulness in the efficient design of decoding simulation 
experiments. 

* This research was supported by the National Aeronautics and Space 
Administration under NASA Grant NSC 5025 at the University of Notre 
Dame in liaison with the NASA Goddard Space Flight Center. 



I. INTRODUCTION 


In this report, we give the results of comparing, by their performance 
on a simulated deep-space channel, virtually all of the rate one-half, 
binary, convolutional codes of rate R = 1/2 and "equivalent constraint 
length" = 24 which have been proposed for use with sequential decoding 
on the deep-space channel. These values of R and K^, are those selected 
for the International Ultraviolet Explorer (lUE) spacecraft. This report 
supersedes an earlier report [1] which gave a similar, but less extensive, 
comparison of = 32 codes 
in the lUE spacecraft. 

Section II contains the definitions of terms used herein and the 
descriptions of the various design options for the convolutional encoder/ 
sequential decoder (CE/SD) system. In Section III, we describe the simulation 
procedures employed and list the various codes which were compared. We then 
give, in Section IV, our recommendations for the code to be used in the lUE 
CE/SD system. We give consideration also in Section IV to many factors which 
should be considered in the development of future coding standards for deep- 
space communications. 

Because of the very small decoding error probability for sequential 
decoding, even at the rather short = 24 constraint length, it was necessary 
to give careful consideration to the statistical significance of the data 
obtained. A novel method to characterize this statistical significance was 
developed and is described in Section III of this report . 


; Kg = 32 having earlier been considered for use 
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II. CODE PAI^IETERS AND DESIGN OPTIONS 


A. Convolutional Coding Terminology 

We give here a brief description of a convolutional encoder for use 
with sequential decoding. Whenever possible, we use the same nomenclature 
as in the preliminary coding standard [2] developed at the Goddard Space 
Flight Center. 

Figure 1 gives tlie general block diagram of a rate R = 1/2, binary, 
convolutional encoder. The encoder is completely specified by its 



Fig. 1. A General R = 1/2 Binary Convolutional Encoder. 


generat or polynomials G^(D) and G^(D), where 

Gj^(D; = + a^^D + ... + ^ 

and 

G^CD) = bp + b^D + ... + 

Eacli binary digit a^ (or b^) is 1 or 0, respectively, according as to wlielher 
or not there is a connection from the corresponding sliif I -regist er stage to 
the corrcsi onding modulo-tvo adder. lor ex.ample, .ij = 1 would signify a 
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connection from the second shift-register stage to the upper modulo-two adder 
in Fig. 1, whereas = 0 would signify the absence of a connection from the 
second shift-register stage to the lower modulo-two adder In Fig. 1. It is 
also customary to speak of the generator functions which are the binary 
K-tuples 


1 ' ®0^1 ‘ ■ ^K- 

2 " *^0^1 • • • '^K-l 


The constraint length , K, is the lengtli of the shift-register in Fig. 1, i.e. 
the number of Information bits in the span which, at any clock instant, 
determine the two encoded bits formed by the modulo-two adders. The infor- 
mation sequence to be encoded is represented by the polynomial 

1(D) = 1q + i^D + ... + 


where i^ is the information bit present in the leftmost shift-register stage 
at clock instant t. Note that, at each clock instant, one Information bit 
enters the encoder shift-register but two encoded bits are formed. These 
encoded bits are multiplexed into one stream to form the encoded sequence 
which is denoted by T(D) in Fig. 1. The quantity L is the frame length in 
information bits . 

Whenever sequential decoding is used, the information sequence is follov; 
ed by a tail of T zeroes so that a total of L + T clock instants are used 
by the encoder to form the entire encoded f rame . This segmentation into 
"frames" is required to permit independent processing of frames so that tlie 
sequential decoder can move on to the processing of the next frame whenever 
it encounters a received frame whlcli causes excessive decoding comi>utal ion ; 
such frames thus result in do 1 et ion of the frame ratlu*r tlmn <i ecodinf’ l^rro rs. 

The choices a^ = b^^ = 1 are always made, so that tliere is always a 
connection from the first r.hi 1 t-register stage to botli moduU)-two adders in 
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Fig. 1. The code is said to be systematic when ••• 

i.e., when the encoded sequence from the upper modulo-two adder in Fig. 1 
is Just the information sequence itself (the upper adder is then, in fact, 
not necessary). Because the information sequence appears directly in 
alternate positions of the encoded sequence, it can be trivially recovered 
from the latter in the absence of errors. In a non-systemat Ic code, such 
simple recovery of the information sequence is not possible. The code is 
said to be quick look-in (QLI) if 62 ( 0 ) = D + G^(D), i.e., if a^ = b^ for 
i ^ 1 but a^ ^ bj^. For a QLI ';ode [3], G^(D) + G 2 (D) = D so that the 
information sequence can be recovered from the encoded sequence, except for 
an unimportant delay of one clock unit, by the simple circuit shown in Fig. 2. 



Fig. 2. Recovery of the Information Sequence from 
Llic Error-Free Encoded Sequence for a QLI 
code . 

Simple recovery of tlie information sequence from the transmitted sequence 
is desirable both for encoder check-out and for extracting tlic engineering 
data from the "hard-dccisioned" received sequence without the necessity of 
full decoding. The probability of error in tlie latter case is the minimum 
possible, viz. tlie error probability in the hard-decisions themselves, for 
systematic a)des. For QLI codes, this error probability is at most twice 
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itmonucmnm of tup 

Olilij’WiVL PAUK is Pixif 


this absolute minimum [3], and is the minimum for nonsystcmatic codes. 

B. Implementation Options and Equivalent Constraint Length 

It is not generally realized that there are several options available 
to the CE/SD system designer as to the choices of K and T even after the 
frame paramete rs liave been frozen . To see this, we sl>ow in Fig. 3 a con- 
ceptual view of the complete encoded frame in an R = 1/2 CE/SD system. 

(We neglect the "synch patterns" which in general are appended to the encoded 
frame to form the entire transmitted frame. The active port? i_n of the encoded 


Active portion 

1 

Passive portion 

2L Bits 

1 P Bits 

1 

1 


Fig. 3. Structure of a Complete Encoded Frame for R = 1/2 

frame consists of the 2L encoded digits formed by the modulo-two adders in 
Fig. 1 during the L clock instants in which information bits enter the shift- 
register. The passive portion of the encoded frame consists of P further 
encoded bite formed during the T clock instants in which the tail of zeroes 
is inserted into the encoder. The total encoded frame length , F, is then 
F = 2L + P. 

We now consider the options still available to tlie CE/SD designer after L 
and P have been frozen. 

1 . Conventional Option 

We first describe the manner in which virtually all past CE/SD systems 
have been implemented, whicli we shall call the "conventional option." Tills 
option is: 

Conventional option: Choose T * K - 1 = P/2. 

In this option, the passive portion of the frame is just the 2T encotlcd bits 
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formed during the T clock instants in which the tall of zeroes enters the 
encoder. This option has the advantageous feature that the encoder may 
Immediately begin encoding the next frame after completing the last because 
when the first information bit of the next frame enters the encoder, it 
automatically finds only O's residing in the remaining T = K - 1 shift 
register stages. 

We now define the equivalent constraint length, K , of any CE/SD option 

to be 

= r/2 + 1. 

Onr rationale for this definition is that 

= K (conventional option), 

and that, sin. c: P completely determines the actual "overhead" in the encoded 
frame needed to effect segmentation of frames, only P should be used as the 
measure of "constraint length" within the frame. 

2. A-Undersized Tail Option 

There is no reason to choose T so large as K - 1 except for the con- 
venience of the "automati': initialization" (as just described for the con- 
ventional option) for encoding of subsequent frames. Tlius, we define: 
A-Undersized tall option: Choc eT=K-l-A= P/2. 

(A, some positive Integer). 

In this option, the passive portion of the frame again consists of the 2'J' 
encoded digits formed as the tall of T zeroes is Inserted into the encoder. 
From tlie definition of Kj, , we see that 

Kj. “ K - A (A-unders i Zed tall option). 

This option lias the minor disadvantage that, In order to begin encoding 
of the next frame, it is lu-cessary to "stuff" zeroes Into tin* A rightmost 
stages of the encoder shift-register immedi.itcly upon conclusion of the 
encoding of the frame. It li.e. , liowever, two m.ijoi advantages over the 
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conventional option. First, the frame dccodlnR error probability , P^, with 
sequential decoding Is in general proportional to the frame length in 
Information bits, L, when K » T + 1 as in the conventional option. However, 
fixing P (and thus also fixing T in the A-undersized tall optlo"), we have 
shown elsewhere [A], [5] that decreases rapidly as K is increased up to 
the point where A = log 2 L where no further improvement results. Moreover, 
for 


A > log^L, 

the resulting P^ is independent of L. We can now sec a second advantage 
in the possibility to standardize encoder design . By choosing one good 
code with a large value of K, we will have a near-optimal CE/SU system over 
a very wide range of telemetry formats, viz., all those such that 
K P/2 + log^E + 1. 


3. Systematic Partial Tail-Suppression Option 

This option can be used only with systematic codes and its possibility 
was first reported by us quite recently [6]. It rests upon tlie obvious 
fact that, since the upper adder in the circuit of Fig. 1 for a systematic 
code emits only zeroes as the tall of T zeroes is inserted into the encoder, 
it is both unnecessary and wasteful to multiplex these known zeroes into 
the encoded sequence. Thus, in this option, the passive portion of the 
encoded fr.nme consists of only the P = T bits from the lower adder in Fig. 1 
during the T clock instants th.it the tall of zeroes is inserted. 

Systematic partial tai ’.-supprer.sion option: Clioose T = K - 1 =* P. 

For this option, we see from the definition of Kj. tliat 
Kj. “ (K - l)/2 ^ 1. 

In this option, we again have T « K - 1 so that automatic initialization 
for successive encoded i ranes is achieved .is in the conventional option. 
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The main advantage over che conventional option, when nonsystcmatic codes 
are used with the latter, is the simpler and mure reliable recovery of the 
information sequence from the hard-decisioned received sequence, although 
this advantage is only a factor of 2 when QLI nonsystcmatic codes are 
considered. There are two minor disadvantages of this option compared to 
the conventional option. First, because only one encoded digit is formed 
at each encoder clock instant during the passive portion of the frame, the 
encoder clock must double its speed when it enters the passive portion of tl>e 
frame. Secondly, the frame decoding error probability, P^, may be slightly 
increased. It has long been known [3] that, for the same K, the best non- 
systematic code gives a far smaller than the best systematic code. At 

R * 1/2, however, the K in this option for the same K (i.o., same P) is 
related r.o the K in the conventional option as 

K - 1 - 2(K -1) 

sys nonsys 

or 

K = 2K 

sys nonsys 

The available evidence indicates that, when their constraint lengths arc 
so related, the best systematic code gives about the same P^ as the best 
nonsystcmatic code [3), but with a sliglit advantage for the nonsystcmatic 
code (7). 

A. Systematic Partial Tail-Suppression with A-undersized Tall Option, 
lliis option is that obtained by combining the two unconventional features 
of the two previous options, viz. removal of known zeroes from the encoded 
sequence and use of a tall shorter than that required to automatically 
initialize the encoder between fr.imes. Tills option cr.n thus be described as: 
Systematic partial tail-suppression 

with A-underslzed tall option: Oioose T“K-1-2A»P 


(A, some positive integer) 
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We Bce for this option that OKIoINaIi PAGli IS PO(ii< 

- (K + l)/2 - A. 

Because this option is a combination of the "unconventional" features 
of the two previous options, it shares both the advantages and disadvantages 
of each. On the minus side, it requires stuffing of zeroes to initialize 
the encoder between frames and also requires the encoder clock rate to double 
in the tail. One would expect also that, for the same and A, the frame 
decoding error probability, P^, would be slightly greater than for option 
2 with a comparably good nonsystematic code; however, this drawback is probably 
not significant. [It should be noted, however, that the constraint lengths 
K in these two options, for the same K^, are related as 

K - 2K - 1. 

sys nonsys 

As will be seen in Subsection C below, this rules out option 4 as a viable 
option because the required exceeds the value of K for which good 

systematic codes e known at the present time .] On the positive side for 
option 4, the resulting would be Ic s than in option 3 and independcnc 
of the frame length in information bits, L. Moreover, one encoder could be 
standardized aid used for all applications where its constraint length K 
satisfies 

P + 2 log.2 h + 1 

Finally, because the encoder is systematic, the information sequence can he 
simply extracted, with maximum reliability, from the hard-decisi oned 
received sequence. 

C. Summary, by Option, of the ■ 24 Codes Compared In This Report. 

The equivalent constraint length, has been specified at - 24 

(or the lUE spacecraft CK/Sl) system, l.c., P - “ 1) * 40 bits have lieen 
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allocated for the passive portion of the encoded frame. 

In the conventional option 1, the specification P ■ ^6 requires that 
a code with K ■ 24 be chosen. 

In the A-undersized tail option, the specification P ■ 40 requires that 
a code with K ■ 24 + A be chosen. V.’c have specified further that A ■ 24 for 
several reasons. First, the resulting K - 48 is about as large a constraint 
length as that for which good nonsystomat ic convolutional codes are presently 
known. Moreover, the corresponding A ■ 24 is large enough to accomodate 
as wide a range of frame longtlis, L, in information bits and allocated 
passive flame length, P, as is likely to be considered in any future deep- 
space applications , so that a standardized K ■ 48 encoder is robust enough 
for all likely future applications. Finally, K ■ 48 is a convenient con- 
straint length for software sequential decoding as it is well-matched to 
most computer wordlengths. Hcr.cc; consider hereafter only K ■ 48 codes 
for option 2. 

In the systematic partial tail-suppression option 3, the specification 
P ■ 46 requires that a systematic code with K ■ 47 be chosen. 

In the systematic partial tai 1 -suppression with A-undersized tail 
option 4, the specification P •• 46 requires that a systematic code with 
K ■ 47 + 2A be chosen. The same considerations as Just discussed for option 
2 suggests that A ■ 24 would be a sensible choice. However, the resulting 
K ■ 95 is far beyond the range for which good systematic codes art presently 
known. In fact, the systematic codes reported earlier [8] with K jc 61 are 
the longest good codes presently known. The K C . code would allow the 
choice A ■ 7 in the IL'L spacecraft CK/SD, but this rather small A would 
accomodate a very smal 1 additional range of fr.ime pariimeters so that standard- 
ization at tills value would be unwise. Because such standardization Is 
perhaps the most attr.icilve fe.iture of option 4 vls-a-vls option 3, 



11 


Its impossibility with the present state of the art in code construction 
has led us to rule out option A as presently unvlablc. 

In Table I, we summarize the types of convolutional codes, by option, 
which we shall consider as candidates for use In the lUE spacecraft CK/SD 
system. 


Option 

No. 

Option 

Name 

1 

K 

Remarks 

1 

i 

Conventional 

1 

2A 

2A 

Viable 

2 

(A“2A)-Undcrsizcd Tail 

2A 

A8 

1 

Viable 

3 

Systematic Partial Tail -Suppression 

2A 

m 

Viable 

A 

Systematic Partial Tail-Suppression 

2A 

95 

Unvlable at 


with (A-2A) Undersized Tail 

, J 


Present 


Tab ‘ i Required constraint length K, for convolutional 
codes with “ 2A , by implementation options. 
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III. SIMULATION RESULTS AND COMPARISONS 

A. Codes Selected for Comparison 

We now list, by option, all of the codes used in the simulations to be 
reported here. These codes comprise all of the known convolutional codes, 
with the required constraint lengths, that appear attractive for use in 
deep-space communications. 

1. Codes for the conventional option. 

a. The Massey-Costello (MC) quick-look- in (QLI) code with K = 24. 
( MCQLI-24 code ). 

This code is obtained by shortening the K = 48 QLI code 
given by Massey and Costello [3]. This K = 24 code is the recommended co(ie 
for this length in the preliminary coding standard [2], 

b. The first Johannesson (J) QLI code having an optimum distance 
profile (ODP) with K = 24. ( lJQLIODP-24 code ). 

This is the first of two K = 24 QLI codes found by 
Johannesson [9] having the optimum distance profile (ODP^ property which 
will be discussed further in Subsection B below. 

c. The second Johannesson QLI code having an ODP with K = 24. 
( 2JQLIODP-24 code ). 

d. The Bahl-Jelinck (BJ) K = 24 complementary code ( BJ-24 code ) . 

Complementary codes were defined by Balil and Jelinek [10] 
as codes for which ~ ~ ^-1 ~ ^K-1 ~ ^ ^i ^ i ^ < K - 1. 

This is the best K = 24 code in this class as found by exhaustive scarcliing 

[ 10 ]. 

e. The K = 24 quadratic residue code. ( QH-24 code ). 

This is a code wliich was found by Massey, Costello raid 



Justesen [11] using a technique for constructing convolutional codes from 


known cyclic block codes. This code actually has K = 23 since = b 2 ^ = 0* 

f. The K *= 24 optimally-truncatable code ( OT-24 code ). 

This code, which was found during the progress of the 
research reported here, will be described more fully in Subsection B below. 

2. Codes for the (4=24)-undersized tail option. 

a. Tlie Massey-Costello QLI code with K = 48. ( MCQLI-48 code ). 

This K = 48 code is the recommended cod', for this length 
in the preliminary coding standard [2]. (However, its use as here at = 24 
in the undersized tail option is not considered in the preliminary standard.) 

b. The Johannesson QLI OOP code with K = 48. ( JQLIODP-48 code ). 

This code was reported in [6] which contains extensions 
of Johannesson 's earlier v^ork [9]. 

3. Codes for the systematic partial tail-suppression option. 

a. The adjoinc [12] of Bussgang's optimal systematic code [12] as 
extended by Lin and Lyne [13], and later Forney [14], to a good, 
but sub-optimum, code at K = 47. ( BLLF-46 code ) . 

This is the K = 48 systematic code recommended in the 
preliminary coding standard [2]. Its use by shortening to K = 47 actually 
gives a code with K = 45 because = b^^ = 0 in the generator. 

b. Johannesson ' s systematic ODP code with K = 47 ( JSODP-47 code ). 

This code was reported in [6] in the extension given there 
of Johannesson' s earlier work. Tliere are no otlier K = 47 systematic codes 
known which appeared promising enough to be considered in this comparison. 

In Table II, we summarize for easy reference the ten Kj. = 24 codes that 
were used in the simulations to be reported here. The generators are listed 
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in octal form in the manner that *= 111* Oil, Oil, 101, Oil, Oil, 110, 111 
(binary, with con aas shown to indicate segmentation for conversion to octal) 
is sho\TO as = 73353367. Also shown in Table II is the free distance , 
dfree* code which is the minimum Hamming distance between two complete 

encoded sequences caused by different information sequences. The minimum 
distance , d 22 » which is the minimum Hamming distance between two encoded 
sequences during the first K = 24 encoder clock instants caused by infer- 
mation sequences that have different values of 1 q, is also given in Table II. 


Code 

No. 

Code 

Name 

G1 

(octal) 

G2 

(octal) 

^^23 

*^free 

la 

MCQLI-24 

73353367 

53353367 

9 

17 

lb 

lJQLIODP-24 

74042417 

54042417 

11 

18 

Ic 

2JQL10DP-24 

74041567 

54041567 

11 

19 

Id 

BJ-24 

51202215 

66575563 

10 

24 

le 

QR-24 

77441232 

54502376 

10 

20 

If 

1 

OT-24 

75105323 

55105323 

10 

19 

2a 

MCQLI-48 

73353367 

67373553 

53353367 

67373553 

9 

>19 

2b 

JQLIODP-48 

74042402 

07121635 

54042402 

07121635 

11 

>_19 

3a 

BLLF-47 

40000000 

00000000 

71547370 

13174650 

9 

19 

3b 

JSODP-47 

40000000 

00000000 

67114545 

75564666 

11 

17 


Table II: 


Description of K = 24 chosen for comparison 
h 

(Where d, is unknouTi , the best known lower 
free 

but these lower bounds are not expected to be 

Instance, it is likely that d, is at least 

■’ free 

and 2b.) 


in simulations, 
bound is given, 
tight. For 
30 for codes 2a 


RKPRODUrmiLITY Oi Vlif.: 

OHl'v'.lXAI. iMOb’ li* [‘Out.’ 
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B, Distance Properties of the Selected Codes 

The distance properties of a convolutional code can be used to make a 
fairly reliable estimate of how that code will perform with sequential de- 
coding. We now consider the distance properties of the above codes with a 
view to correlating their performance on the deep-space channel with these 
properties. 

The i-th order column distance , d^ (i = 0,1,2,...) is defined to be the 
minimum Hamming distance between two encoded sequences over the first i + 1 
time instants [2(i+l) digits for rate 1/2 codes] resulting from information 
sequences with different values of i^. Because of the linearity of the code, 
d^ is equal to the minimum Hamming weight of an encoded sequence over the 
first i + 1 time instants resulting from an information sequence with 


ip = 1. The free distance , may be expressed as 

d^ = lim d, 
i-Ho i 

and is well-known to be the maiii determiner of error probability where 
sequential decoding is used [3]. The distance profile , [d^ ,d ,d 2 , • . . d^^^j , 

has been shown [9] to be the property of interest as far as computational 


performance with sequential decoding is concerned. The distance profile ^ 

is said to be superior to the distance profile if d^ > dj^ for the smallest 

i such that d^ =f d!^; we denote this superiority by writing ^ . We define 

the distance index i^ of a code to be the smallest integer i such that 

dj = i. Because d„ = 2 for all R = 1/2 codes having a„ = b^ = 1 , i» is the 

first distiince index of interest. Because d. < d^ for all i, i, is 

1 — free d. 

free 

the largest distance index of interest. The distance indices ij for 
1 = 2, 3,..., dr of course uniquelv determine the distance profile d, but 

j > > » f rce 

in general they also provide further information, viz. that about tlie column 


distances d^ for i > M. 
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In Table III, we show the distance indices for all of the codes to be 
compared in our simulations of K « 24 codes. For codes Id, 2a and 2b, the 
distance indices are shown wherever we have found their value, but the time 
requirements for computing higher distance indices were so great as to pro- 
hibit generation of th'J full set of distance indices. All of the distance 
indices shown were obtained by using a sequential decoding program, as was 
first suggested by Forney [14], which was programmed in assembly language 
for the IBM 370/158 computer. In this procedure, for a given threshold 
T = j-1, one explores all encoded sequences beginning with i^ = 1 until 
their Hamming weight exceeds T. If the longest path of weight T = j-1 is 
i + 1 branches in length, then d^ = i. The program also checks the encoder 
state after each extension of an encoded path whose Hamming weight is T = j-1. 
If the all-zero state is ever encountered, then = j-1 and the program 

is stopped. Otherwise, the threshold T is increased by 1 and the search 
continued until either a zero state is reached or the predetermined limit 
on CPU time is exceeded. The values of d^ given in Table II were all 
found in this manner, except that for the BJ-24 code whose was found 

in [10]. The CPU time required for each code in Table II with d^ >19 
was approximately 15 minutes. 

It is interesting to compare the codes in Table II in terms of their 
distance properties. In terms of their free distance, the ordering by 
quality of the codes is 

(2b) I (2a) I (Id) > (le) > (Ic) = (3a) = (If) > (3b) = (la) 
[ordering by free distance] 

where the codes are sliowi in order of decreasing quality. In terms of their 
distance profile (whore we treat all codes as having M = 47, the memory of 
the code with greatest memory), the ordering by quality of the codes is 
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Code No. ! la lb Ic Id le If 2a 2b 3a 3b 



Table III: Distance properties of the codes chosen for comparison in simulations, 

(ij is tl)c smallest i such that = j.) 
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(2b) - (3b) > (Ic) > (lb) > (3a) > (Id) > (If) > (la) > (2a) > (le) 

[ordering by distance profile]. 

As we shall see below, these two orderings substantially agree with the 
ordering in performance with respect to decoding error probability and 
computation, respectively , for these codes on the simulated deep-space channel. 
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C. Description of the Simulation 

The channel chosen for simulation was a white additive Gaussian noise 

channel (l.c., a deep-space channel) with an energy per transmitted bit, 

E, to one-sided noise power spectral density, N , ratio of E/N = 1 or 

o o 

0 db. Because the code rate is 1/2, this corresponds to an energy per 
information bit, E, , to N ratio of E, /N = 2 or 3 db. The received digits 
were quantized to 3 bits using the quantization scheme suggested by Jacobs 
[15]. The result of this quantization is to convert the Gaussian channel 
to a 2-input, 8-output discrete memoryless channel (DMC) whose transition 
probabilities are as given in Table IV. This DMC is tlie actual channel 
that was simulated, the various transition probabilities being obtained by 


-^--^output 

input 

0 

O' 

O' ' 

O' ' ' 

1'" 

1" 

1' 

1 

0 

.434 

.197 

.167 

.111 

.058 

.023 

.008 

.002 

1 

.002 

.008 

.023 

.058 

.111 

.167 

.197 

.4 34 


Table IV : Transition probabilities for the DMC obtained by 3-bit quanti- 

zation of the deep-space channel with E, /N = 3 db for R = 1/2 coding. 

b o 

use of a standard random number routine. The cut-off rate, R , of the 

o 

Gaussian channel (whicli is sometimes also denoted as R ) was .5481 , 

comp 

whereas that of tlie resulting DMC is .4913 so that the loss due to quanti- 
zation is 0.48 db. Note tliat R^ = .4913 is slightly less tlian tlie code rale 
R = .5000, which indicates that this DMC is actually somcwliat noiser than 
the channels £0 wlilch the CE/SD sys tern would actual ly be used . This "noisier 
than usual" situation was chosen for simulation so tliat the probability of 
a decoding error would be large enough so that the number of decoded frames 
required to determine the decoding error probability would correspond to 
practical amounts of computing tiitic. 
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The Fano sequential decoding algorithm [16] was used for decoding. 
The metrics used were as given in Table V. 







mm 

B 

B 

IB 

0 

B 

B 

3 

B 

-8 

-20 

-34 

-59 

1 


-34 


-8 

0 

3 

4 

4 


Table V ; Metrics used for sequential decoding for the DMC given in Table IV. 

The threshold increment A used in the Fano algorithm was 32. A frame size 
of L = 256 information bits was selected. A frame was declared to be 
erased when 100,000 computations were performed without completion of the 
decoding. This erasure limit of 390 computations per decoded information 
bit is well above what one would generally choose in practice, the reason 
for this choice being again to ensure a significant number of frames with 
decoding errors (rather than merely erasing most of those frames.) The 
number of computations used for an erroneously decoded frame was recorded 
so that one can determine how would change if a smaller erasure limit 
on computation were chosen. 

For all codes selected for comparison, 10,000 frames were decoded in 
the simulation. The same 10,000 received frames were decoded for each code 
so that any differences in computation and/or error probability w’ould be 


due to the codes themselves. 
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Code No. 


Code Name 


la 


«a- 

CM 

I 

M 

O' 

u 

T. 


lb 


CM 

I 

CL, 

a 

o 

hJ 

o 


Ic 


CM 

I 

ou 

Q 

O 

M 

O' 

(N 


Id 


le 


If 


CM 

I 

—i 

cq 


CM 

I 

a: 

o 


>T 

CM 

H 

O 


2a 


00 

I 

M 

►J 

O' 

u 

X 


2b 


00 

I 

CL. 

Q 

O 

M 

_) 

O' 


3a 


C'- 

«J 

I 

Cl. 

hJ 

ca 


3b 


r.. 

I 

CL, 

Q 

O 

(,) 


ilfFrarae errors 
//Bit errors 
//Erased frames! 


A 

A8 

35 


A 

59 

2A 


2 

36 

23 


0 

0 

25 


0 

0 

27 


5 

116 

25 


1 

38 

26 


0 

0 

22 


1 

13 

17 


0 

0 

2A 


//Frames with: 


C £ AOO 

1058 

1122 

1118 

1103 

C 500 

3029 

3283 

3A2A* 

3352* 

C _< 650 

51A5 

5575 

5552* 

5A26* 

C £ 850 

6681 

7029 

7012 

6965 

C < 1,200 

7852 

8136 

8115* 

8051* 

C < 1,700 

8558 

875A 

875A* 

8696* 

C _< 2, AOO 

8992 

9109 

9119* 

907A* 

C 3,500 

9298 

9385 

9397* 

9363* 

C _< 5,000 

9A97 

9563 

9550 

95A0 

C ^ 7,000 

9637 

9671 

9672* 

9652* 

C ^ 10,000 

9737 

9762 

9765 

97A1 

C _< 15,000 

9823 

9826 

9835* 

9828* 

C _< 20,000 

9859 

9873 

9872 

9871 

C 30,000 

990A 

9916 

9915 

9913 

C _< 50,000 

99A1 

9950 

99A8 

99 A 8 

C ^ 100,000 

9965 

9976 

9977 

9975 

C in fr.ames 

36 , 369 

1 1.1 R5 

70,213 


with errors 

12,A91 

16,3A6 

9,7A0 



933 

23,A10 




A8,065 

85,658 




1037 

1100 

1063 

1123 

611 

602 

2918 

319A 

2987 

3288 

2653 

2689 

A961 

53AA 

5175 

5588 

5135 

5268 

6518 

6838 

6695 

7013 

682A 

68A9 

7755 

7957 

7860 

81A2 

800A 

8080 

8507 

8627 

8561 

8756 

8702 

87AA 

8971 

90A7 

8983 

912A 

9098 

9106 

9300 

93A6 

9290 

9A02 

936A 

9382 

9507 

953A 

9A85 

9552 

9553 

95A6 

9633 

9665 

9613 

9678 

9681 

9672 

97A5 

9756 

9729 

9771 

97 7A 

9773 

9827 

9831 

9822 

9833 

9855 

9831 

986A 

9876 

9869 

9870 

9881 

9867 

9897 

9913 

9902 

9922 

9919 

9920 

99A0 

99A6 

99A7 

99 A 8 

99A2 

9950 

9973 

9975 

997A 

9978 

9983 

9976 


73A 

A9,AA9 1 

1 

17,922 1 


62,585 

37,]A9 

2,969 

3,0A7 


Table VI: Results of Decodinp 10,000 frames of lenj;lh I. * 256 information bits 

on the simulated dl•el)-^ipace channel wltli an K, /N of 3 db. (V.ilucs 

b o 

marked * were interpolated from nearest d.ita points obtained in tlie 
simulation on the .assumption that the computation is I’.are to-d 1 t r 1 butrd • ) 



















22 


D. Simulation Results 

In Table VI, we show the results obtained from decoding the same 10,000 
frames of length L ■ 256 on the simulated deep-space channel with an 
of 3 db. Because no decoding frame errors were made with four of the codes, 
additional frames were decoded for these codes with the results shown in 
Table VII. 


Code No 

Id 

le 

2b 

3b 

Code Name 



00 





1 





Pu 





Q 

1 

i 

-a- 


O 

Pu 


CN 

<N| 

M 

Q 


1 

1 


O 




o 

to 


PQ 

O' 



(/Additional 





frames 

40,000 

40,000 

40,000 

10,000 

decoded 





//Frame 

0 

2 

0 

2 

errors 


1 

1 


//Bit errors 

0 

20 

0 

41 

C in frames 


23,864 

1 

9,249 

with errors 


45,786 

- -1 

57,636 


Table VII ; Results of additional decoding simulations for those •'.odes 

wliich gave no d oding errors on the first 10,000 decoded frames. 

Approximately 15 minutes of CPU time on the IBM 370/158 system were required 
for each 10,000 decoded frames. .Neither code Id (the BJ-24 code) nor code 
2b (the JQLIODP-^8 code) yielded any decoding errors on a total of 50,000 
decoded frames. It was not feasible to decode enough additional 1 rames to 
distinguis)) wliich code .ictually gives the better decoding error prob.ihi 1 i ty . 

As regard.s decoding error probability as determined by simulation, it 
can be seen from Tables VI and VII that the codes rank In quality as: 
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(2b) - (Id) > (le) > (3b) > (3a) > (2a) > (Ic) > (la) > (lb) > (If) 
where we have resolved ties In the number of frames decoded In error by 
considering the number of Information bits decoded In error. We note that 
this empirical ranking Is substantially the same as the ranking by 
given In subsection B above. 

In comparing the frame decoding error probability, , one should also 
give attention to the number of computations done in the erroneously decoded 
frames as this can be used to determine the change in should a different 
frame erasure limit be used. For instance, we see from Table VI that if 
the (more realistic) limit of 10,000 computations for a decoded frame is 
used, codes la, lb and If would yield 1, 1 and 3 erroneously decoded francs, 
respectively, compared to 4, 4 and 5, respectively, for the 100,000 
computation limit. 

As regards computational performance as determined by simulation, it 
can be seen from Table VI that the codes rank in quality as 

(lb) > (2b) > (Ic) > (Id) > (3b) > (If) > (3a) > (2a) > (la) > (le) 
where we have, somewhat arbitrarily taken the number of frames decoded with 
850 or fewer computations as the Indicator of quality. We note that this 
ranking is in substantial agreement with the ranking by distance profile, 
d, given in subsection B above. In particular, we see the excellent compu- 
tational performance of all of the optimum distance profile (OUP) codes. 

A word should be said about the apparently poor computational performance 
of the two systematic codes (codes 3a and 3b) when the computation is sm.ill. 
Tlic reason for this effect is that a tail of T - 47 branches (with only one 
digit per branch) is included in the frame for these codes, compared to only 
T ■ 23 branches (with two digits per branch) for all the other codes. 


2A 


Because "computation" is measured in terms of branches searched, the T • 47 
systematic codes necessarily require a small number of additional computa- 
tions per frame, compared to the T 23 nonsystcmatic codes, to account for 
the computation done on the extra tail branches. However, as the computa- 
tion per frame increases, the additional computation becomes a negligible 
part of the total computation (which is the reason we compared computational 
performance at 850 computations or less per frame rather than some smaller 
number. ) 
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E. Jtatlstlcal Slcnlflcance of the Decoding Error Probability aa 
Determined from the Simulations. 

The very small number of frame decoding errors, between 0 and 5 inclusive 
for the best and the worst codes, makes it essential to consider the statis- 
tical significance of the frame error probabilities determined thorc''rom. 

We now present a new method for testing this significance, a method which 
it appears could be of rather general usefulness in decoding simulation 
studies. 

Suppose that N frames are decoded for a particular code and that the 

true frame decoding error probability is P^. Because decoding of f rames is 

done independently , and because N >> 1 while << 1, the number of frame 

errors, X, is a Poisson-distributed random variable whose mean and variance 

are both equal to 

X - NP . 
e 

Let X be the observed value of X for the particular N frames actually decculecl. 

Tlicn X is the best estimate of both E(X) *• X and Var(X) ■ X, while x/N is 

the best estimate of P ■ E(X)/N. The problem is that, for >: small (sav, 

e 

X < 10), the estimate of Var(X) is quite inaccurate so that the standard 

statistical "confidence level" techniques i annot be appli»“d with any validity. 

We now show how this dilemma can be resolved. 

We define the 10U3a con i i dence i nt e rv. tl , ('j. • £fLT f be 

observation x of X, in the m.inner that X and > are the large.-, t .md ^i^.a’ll^'.t 

L L 

numbers, respectively, sucli that 


P(X < 

x) P, 

all X < A j , 

PCX > 

X 
1 V 

all X _> X^j. 
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This Is equivalent to saying that, had the code possessed an actual \ smaller 
than the probability of observing fewer than the number of x of frame 
errors actually observed would have been at least 6; and that, had the 
actual X been larger than X^, the probability of observing more than the 
number of x of frame errors actually observed would have been at least B. 

Because X is Poisson with mean and variance equal to X, the frequency 
distribution for X, i.e., the probability of the event that x equals i, 
is given by 




Now, letting F^(X) be the probability that X ^ j. 


F/X) 

J i=0 


I X^ -X 


we have 


It now follows, from the definitions of X and X , that 

Lj u 


and 


Vi<V - « 


F^(X„) = 1 - S. 


For a given x, these equations may be readily s'-lved, e.g. , by the Newton- 


Raphson method, to determine and In Table VIII, we give the 90% 

confidence intervals, so determined, for 0 < x < 5. 


X 



0 

0 

2.30 

1 

.105 

3.83 

2 

.532 

5.29 

3 

1.107 

6.65 

4 

1.76 

7.99 

5 

2.52 

9.06 


Table VIII ; 90% confidence intervals for X given the observed value x 
for a Poisson random variable X. 


■ ■ •; or Tliti 
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Given the lOOSZ confidence interval (X^^, X^) for X, the corresponding 

1000% confidence interval ^e " E(X)/N - X/N is just 

(X. /N, X,,/N). In Table IX, we give the 90% confidence levels for P as 

determined from the simulation results given in Tables VI and VII. 

When the confidence interval for one code lies entirely to the left of 

that for another, then with at least 1000%' confidence we can assert that 

the actual P „ of the former code is smaller than that of the latter. In 
eU 

Figure 4, we have plotted the 90% confidence intervals on a logarithmic scale 
for all of the codes tested. We see from this figure that with at least 90% 


Code No. 

// frames 
decoded, N 

// frame 
errors, x 

! best estimate 

i of P , x/N 
e 

1 90% confidence interval 

! "'eL- ‘’eU> f”-- ’’e 



1 


i _5 _5 

la 

10,000 

! 4 

40 X 10 " 

(18 X 10 80 X 10 




-5 

-5 -5. 

lb 

10,000 

4 

40 X 10 

(18 X 10 80 X 10 

1 

Ic 

10,000 

2 

20 X 10“^ 

(5.3 X 10"^, 53 X lO"^) 

Id 

50,000 

0 

0 

(0 , 4.6 X 10"^) 

le 

50,000 

2 

4 X 10“^ 

(1.1 X 10"^, 11 X 10"^) 




-5 

. -5 -5 

If 

10,000 

5 

50 X 10 ^ ; 

(25 X 10 , 91 X 10 




1 

-5 

-5 „ -5x 

2a 

10,000 

1 

10 X 10 

(1.1 X 10 , 39 X 10 ) 




0 ! 

c 

2b 

50,000 

0 


(0 , 4.6 X 10"^) 




-5 

^-5 „ ..-5. 

3a 

10,000 

1 

10 X 10 

(1.1 X 10 , 39 X 10 ) 




-5 

^-5 -5. 

3b 

20,000 

2 

10 X 10 ^ 

(2.7 X 10 27 X 10 


Table IX: 90% confidence intervals for P as determined from the 

e 

decoding simulations reported in Tables VI and VII. 
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»■ 


3a 

2a 


Id 

2b 


le 


3b 


Ic 


M- 


la 

lb 


If 


10 


10 ' 10 ■ 10 
Fig. 4. Graphical portrayal of the 90% confidence inter'^als for P 


-3 


for all codes tested. (Codes shown by code number), 


confidence we can assert that the frame error probability for codes Id and 
2b (neither of which yielded any errors in 50,000 decoded frames) is truly 
smaller than that of all the other codes, except codes le, 2a, 3a and 3b where 
our confidence in the superiority of Id and 2b is somewhat less. It is a 
remarkable property of Poisson-distributed random variables that such a 
small number of error events (at most 5 for any code tested) can be so highly 
significant statistically. 

It probably should be pointed out that, although 256 information bits are 
decoded in each frame so that there are 256 times as many bit decoding 
decisions as frame decoding decisions, one cannot assert greater statistical 
confidence in the observed decoding bit error probability than in the observed 
decoding frame error probability. The reason of course is that the decodings 
of bits within a frame are highly dependent so that one has no more independent 
bit decoding decisions from which to infer probabilities than one has inde- 


pendent frame decoding decisions. 


IV. DISCUSSION AND RECOMMENDATIONS 


A. Performance Summary for Each Code 

1. Codes for the conventional option, 
a. The MCQLI-24 code, (Code la). 

Although this is the recommended ! * 24 in the preliminary coding 
standard, its frame error probability was tied for worst with lJQLIODP-24 
code among all codes tested, and significantly so (cf. Table IX and Fig. 4). 
Moreover, its computational performance was superior only to the QR-24 code 
(cf. Table VI). Perhaps, however, this is a good place to make the point 
that only codes known to be good were tested . One might say that this 
code is the worst of a good lot! Moreover, it does possess the desirable 
quick-look-in property which is worth some sacrifice in performance (but is 
also a feature of the distinctly better 2JQLI0DP-24 code.) 

Because this code is the K = 24 truncation of the MCQLI K = 48 code 

whose K = 32 truncation is an excellent code at that length [1] , its poor 

K = 24 performance was somewhat unexpected. To explain this phenomenon, 

the free distance for all the distinct truncations of the MCQLI K = 48 code 

with K £ 28 were computed with the results given in Table X. The value of 

d, = 23 for the K = 32 code is also sho\/n. From this table, one sees 
free 

that there is indeed a "soft spot" around K = 24 in the truncations. In fact, 
the K = 22 code is quite likely superior in to the K = 24 code because 
the same free distance is achieved within a smaller constraint length! 


K 

2 

3 5 

6 8 9 

10 

12 14 15 

17 

18 

19 

20 22 23 

24 

25 

26 

28 

32 

^frce 

3 

5 6 

7 9 7 

9 

11 11 11 

13 

11 

13 

15 17 15 

17 

18 

18 

^21 

23 

Table 

A- 

Free 

distance 

for 

truncation 

to 

constraint 

length K, 

for K » 

32 




and all K £ 28 giving distinct truncations, of the MCQLI-48 code. 
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Beyond K > 2V, the codes rapidly Improve; the K - 32 code being already 
quite good. From the data in Table X, we conclude that would probably 
be unwise to recommend truncations of the MCQLI-48 code with K < 32 In 
future coding standards . 

b. The lJQLIODP-24 (Code lb). 

This code gave the best computational performance of all codes tested 

(cf. Table VI). However, its superiority in this respect over the 

structurally-similar 2JQLIODP-24 code was very slight, but its inferiority 

to the latter in P (cf. Table IX and Fig. 4) is substantial. It thus 
e 

appears unwise to recommend this code for any applications . 

c. The 2JQLI0DP-24 code, (Code Ic) . 

The three nonsystematic optimum distance profile codes tested (codes lb, 

Ic and 2b) all gave virtually the same computational performance and were 
superior in this respect to all other codes tested (cf. Table VI), often 
quite significantly. The 2JQLIODP-24 code was also the best of the K = 24 
codes with the quick-look-in feature (codes la, lb, Ic and If) with respect 
to P^ (cf. Table IX and Fig. 4), and significantly so. This code is 
definitely the best of the K = 24 QLI codes studied . 

d. The BJ-24 code, (Code Id). 

This code gave outstanding P^ performance (cf. Table IX and Fig. 4), as 
expected from its large of 24. Its computational performance uas 

slightly inferior to that of the ODP codes (cf. Table VI), but not signifi- 
cantly so. This code, unfortunately, does not have the quick-look-in property; 
but we have shown elsewhere [17] that it can be encoded almost as simply as 
a QLI code although, of course, the extraction of data by an encoder inverse 
from the hard-dccisioned received sequences is much less reliable than for 
a QLI code. However, J_f one docs not insis t on the QLI property , this code 
is clearly the best of al 1 known K ■ 24 codes . 
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e. The QR-24 code, (Code le). 

As we noted earlier, this Is actually a K - 23 code as both of Its K - 24 
generators end In a 0. This code Is primarily Interesting In that It Is the 
only one of all the codes tested which was derived by an algebraic construction 
[11] as opposed to being found by a heurlstlcally-gulded computer search. 
Although Its computational performance Is distinctly Inferior to the ODP 
codes and the BJ-24 code (cf. Table VI), Its performance was surprisingly 
good (cf. Table IX and Fig. 4). However, because It also lacks the QLI 
property (while being inferior in both P^ and computation) as does the BJ-24 
code, this code should not be recommended for any K = 24 application . Its 
superiority in P^ over most of the codes tested does suggest, however, that 
there may well be merit in Investing further effort toward finding algebraic 
approaches to convolutional code construction , an almost totally-neglected 
research area. 

f. The OT-24 code, (Code If) 

The rather uneven quality of the various truncations of the MCQLI-48 
code (cf. Table X) suggested to us the advisability of devising a code 
which would truncate well everywhere. This led us to define the optlmally- 
truncatable OT-K code as the QLI code obtained by the following algorithm: 

Step O i Set 3 q = 1, bg “ 1 , aj^ = 1 and b^^ «• 0. Set 1 * 2, d « 3 

and 1 . = 1. 
a 

Step 2 ; Set a^ = b^ “ 1 and compute the free distance d' of this code 

and the distance index i,,. 

a 

Step 3 : If d' > d or if d' = d and 1^, < 1^, replace d by d’ and go 

to step 4. Otherwise, set 


Step 4 t If 1 “ K-1, stop. Otherwise, increase 1 by 1 and go to step 2. 
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It should be clear that, for K* < K, the OT-K' code is the K' truncation of 
the OT-K code, so that these QLl codes are "nested" In the sane sense as the 
MCQLl codes. 

In Table XI, we give the free distance for all the distinct OT-K codes 
with K ^ 24. Comparison of Table XI with Table X shows indeed that the OT-K 
QLI codes are indeed superior to the MCQLI-K codes, but significantly so 
only for rather small K. However, if one desires a QLI code , in future 


K 

2 

3 

4 

6 

9 

13 

15 

17 

18 

20 

23 

24 


3 

5 

6 

8 

10 

11 

14 

14 

16 

17 

18 

19 


Table XI ; Free distance for distinct truncations of tne K = 24 optimally- 
truncatablc code. 

coding standards , which truncates well for K < 24, the OT-24 code would be 
£ good choice . 

The penalty for both truncating well at all smaller constraint lengths 

and having the QLI property, however, is that one cannot do as well as an 

unconstrained code at some fixed large K. For instance, the OT-24 code 

has d- = 19 compared to d^ = 24 for the BJ-24 code, 
free free 

From Table III, we see that the OT-24 code is superior to the MCQLI-24 
code in both distance profile and free distance. This explains the slightly 
better computational performance of the OT-24 code (cf . Table VI) , but 
raises questions about the observed superiority in for the MCQLI-24 code 
(cf. Table IX); however, as Figure 4 clearly shows, the observed difference 
in for these two codes (la and If) is not statistically significant. 

2. 4-undersized tail option. 

a. The MCQLl-48 code, (Code 2a). 

This code, which is the K ■ 48 code recommended in the preliminary 
coding standard [2], performed (as expected) slightly bettor computationally 
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than its K ■ 24 truncation, code la (cf. Table VI). While its P performnnce 

e 

is certainly good (cf. Table IX), it is somewhat surprising that even one 
frame was incorrectly decoded since the theoretical bounding argument [4] 
suggests that, because of the extra memory, should be reduced by a factor 
of at least “ 10.7 over its K = 24 truncation for which only 4 frame 

decoding errors were observed. The explanation is provided by Table 
the actual error probability of this K = 48 code, for which only x = 1 frame 
decoding error was observed, has a significant probability of being smaller, 
by at least the theoretically-expected order of magnitude, than that 
observed in the simulation. However, because this code is still inferior 
in all respects to the JQLIODP-48 code to be described next, it would be wise 
to rep] ace the MCQLI-48 code in the standard by the JQLIODP-48 code if the 
disruption would not be too great. 

b. The JQLIODP-48 code, (Code 2b). 

This code performed superbly in both computation (cf. Table VI) and in 

error probability (cf. Table IX and Fig. 4). The excellent computational 

performance, which was the best of all codes tested, results of course from 

its optimum distance profile over the full K = 48 branches. The extremely 

low observed P , which tied with that of the BJ-24 code for the best of the 
e 

codes tested, is explained by the expected reduction by at least a factor 
of L/Xj. “ 10.7 over the of its K = 24 truncation because of the additional 
memory of the K = 48 code. 

Tliis is probably the appropriate place to point out that £ truncation 
of an ODP code is al so again an ODP code. Thus, all truncations of the 
JQLIODP-48 code would perform well computationally compared to otlier codes 
at their constraint length. The JQLIODP-48 code would be a good choice for 


the truncatable K » 48 quick-look- in code in future coding standards. 
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3. Systematic partial tail-suppression option. 

a. The BLLF-47 code, (Code 3a). 

This code, which is the K = 47 truncation of the K ■ 48 systematic 

code given in the preliminary coding standard [2], had nearly as good a 

computational performance as the JSODP-47 code (cf. Table VI) to which it 

most closely compares, and is also nearly as good in error probability 

(cf. Table IX and Fig. 4). In fact, as can be seen from Table III, the 

distance profile of this code is only slightly inferior to the ODP code 

and, moreover, this code has a larger free distance than the ODP code 

(19 vs 17). [The paradoxically slightly better P^ of the ODP code derives 

most likely from the fact that it reaches its free distance much more quickly 

(ij^y “ 46 compared to i^^^ = 53 and i^^g = 59 for the BLLF-47 code).] The 

good P^ performance of both codes confirms the conjecture, noted in Section 

II. B. 3 above, that a good systematic code in option 3 performs about as v;ell 

in P as does a good nonsystematic co.'e in option 1, for the same K_. 
e t. 

The "systematic" property is a highly desirable one for many purposes 
(such as the absolutely minimal error probability in recovering data by an 
encoder inverse from the hard-decisioned received sequences, and the absolutely 
simplest "encoder inverse".) The disadvantages in option 3 compared to 
option 1 are the necessity to double encoder speed in the tail, and the small 
amount of extra computation required because of the extra branches in the 
tail. 

The rather slight inferiority of the BLLF-47 code to the JSODP-47 code 
suggests that titere would be only si ight advantage in replacing this code 
by the better JSODP-47 code in future coding standards . 

b. The JSODP-47 code, (Code 3b). 

As Just discussed, this is the best available systematic code with 
respect to both computation and decoding error probability. Its ODP structure 
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guarantees that Its truncations will also be ODP, and hence good. If the 
changeover from the BLLF-A7 code would not be too disruptive, it would be 
the best available choice for the K ■ 48 truncatable systematic code in 
future coding standards. 
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B. ReconunendaClons for the ■ 24 code 


REPRODUCIBILOT OP THE 
ORIGINS PAGE IS POOP 


In the performance summaries for the codes tested, we have already 
made a number of recommendations regarding the selection of a - 24 code 
for the deep-space channel and, more generally, for the selection of rate 
R ■ 1/2 convolutional codes for the deep-space channel. We now particularize 
our recommendations to the choice of a ■ 24 code for the lUE spacecraft. 

1. The best choice for the ■ 24 code would be the JQLIODP-48 code 
(code 2b) used in the (A ■ 24)-undersized tail option (option 2). 

We favor this choice because (a) this code has the best computatioi.al 
behavior of all the codes tested, and a decoding error probability tied with 
the BJ-24 for the best, (b) this same code and encoder could be advantageously 
used at “ 32, or any ^ 48, in later deep-space missions so that its 
adoption now v»ould be a healthy move toward encoder standardization, (c) the 
lone drawback of this code (or any code used in option 2) , namely that the 
encoding shift-register must have all its contents set to zero upon completion 
of encoding of the frame, appears to be very slight, and (4) this code has 
the desirable quick-look-in (QLI) property. 

2. The next best choice would be the BJ-24 code (code Id), provided 
the QLI feature is not considered essential. 

If the conventional option (option 1) is chosen for the lUE convolutional 
coding system, this is definitely the best code to use. Its decoding error 
probability was equal to the best of the other codes tested, and its compu- 
tational performance was nearly as good as any other. The only drawback for 
this code is tiiat it does not have the quick-look-in feature. 

3. If the QLI feature is considered essential, the next best choice 
(after the JQLlODP-48 code) would bo the 2JQL10DP-24 code (code Ic) or the 

K ■ 24 truncation of the JQLIODP-48 code, although the MCQLI-24 code (code la) 


would be acceptable. 
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If both the conventional option (option 1) and the QLI feature are 
chosen for the lUE convolutional coding system, Che 2JQLIODP-24 code offers 
the twin advantages of near-optimal computational performance and near- 
minimal decoding eir r probability. Although the MCQLI-24 code was somewhat 
xnferlor In both of these respects, the fact that various spacecraft have 
already used truncations of the MCQLI-48 code, so that proven encoder hard- 
ware Is presently available, may make this code the practical choice for the 
lUE spacecraft. V.owever, If the JQLI0DP-A8 code (code 2b) displaces the 
MCQLI-A8 code (code 2a) as the truncatable QLI code In future coding standards 
(In keeping with our recommendation In Section IV.A.2.b), then the K ■ 24 
truncation of this code should be used In preference to the MCQLl-24 code. 


POSTSCRIPT 

Following just upon completion of this report, the author received an 
advanced copy of the paper "Further Results on Binary Convolutional Codes 
with an Optimum Distance Profile" by R. Johannesson and E. Paaske to be 
presented at the IEEE International Symposium on Information Theory In 
Ronneby, Sweden, June 21-24, 1976. This paper reports their newly discovered 
K ■ 24 code, namely 

G1 - 55346125 (octal) 

G2 - 75744143 (octal), 

which has both an optimum distance profile and a free distance d^ ■ 25 
exceeding that of any previously known K ■ 24 code (Including the BJ-24 code 
which has d ■ 24.) This code is thus undoubtedly superior both in computation 
and in error probability to the BJ-24, and should be recommended in place of 
the latter code wherever the latter code is recommended in this report. 
Simulations to confirm this superiority are now in progress and the results 
will be reported as soon as they become available. 
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